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COMPRESSIBILITY  IN  STRESS  ANALYSIS* 

By  N.C.  Huang,  E.H.  Lee  and  T.G.  Rogeru 

'  Stanford  University 
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Abstract 

'Almost  nil  the  solutions  in  the  liteianturc  of  visco¬ 
elastic  stress  analysis  have  dealt  with  materla.Ic  which 
exhibit  viscoelastic  behavior  in  shear  but  are  incompress¬ 
ible  or  elastic  in  dilatation.  This  has  been  influenced 
iri  part  by  the  fact  that  accurate  measurement  of  visco¬ 
elastic  characteristics  in  dilaf.ation  are  dlfPlcult,  to 
mak(?,  and  only  now  arc  preliminary  j’csultc  becoming  avail¬ 
able.  Since  inor<'  lirocico  measurement, s  of  dllatatlonal 
I’or.ponsc  can  ho  ■•xpect.ed,  certain  solutions  of  viscoelastic 
stress  analysis  problems  are  presented  to  illustrate  the 
cotivcnicnc''  of  utilising  nioasured  i'olaxation  moduli  in 
dilat,ation  .as  W('‘il  as  siioa!'  in  ova  I  ua  t  In;",  stress  distri¬ 
butions.  Example.;  for  cylindrical  boundaries  arc  j;;lven. 

Thes.o  .;o  I'm,  ions  ai'o'  also  us'M'ul  in  sugsa 'St,  i  up;  '.'xpor  linont  a  I 
configurat  k.ns  fr’oin  wliicli  tin;  r'.-Jaxatjon  modulus  in  dilata- 
t.ion  can  l<;  Vh'ducod  as.  t  is  mia  t  lea  1 1  y  while  avoiding  1‘.h<.' 
d  1  f  f  leu  !  i.  le.-  o  f  dii-ec'.  iiioasuromtni;,  by  su  rround  iti;",  :i  specimen 
wlt.h  fluid  und','!'  pri-’.s'.ure .  M-.'t.aJ  casuii'.s  can  be  i  ncoi'pora  to'd 
_ 

The  results  communicated  In  this  I'opoi’t  wena;'  obtained 
in  Ihe  course  of  .an  invest  i ga L  i  on  conducted  und'O'  Contract 
Nenr'-22',  (on)  by  £  aiii'ord  Uriivenslty  wltli  ‘.no  Office  of 
Naval  R'..'.;e.arch  ,  Washington,  D.  C.  Reproduct. ion  in  whole  or 
in  part  is  [jermilted  for  any  purpose  of  the  United  Si  nt.cs 
Uoverninent  .  Tills  work  i t  o  i.e  in’osont  ed  at,  the  Fourth 
Intern- tiortil  Corisjoss  on  Ra.' ‘ol  0:",v ,  Brown  University, 

Prov  id'-'.nc' ■ .  R.  I..  Au.  ;ust  I'"-.:. 


ON  THE  INFLUENCE  OF  VISCOELASTIC 
COMPRESSIBILITY  IN  STRESS  ANALYSIS* 

V  By  N.C.  Huang,  E.H.  Lee  and  T.G.  Rogers 

\  Stanford  University 

Abstr^t 

^Almost  all  the  solutions  in  the  literature  of  visco¬ 
elastic  stress  analysis  have  dealt  with  materials  which 
exhibit  viscoelastic  behavior  in  shear  but  are  incompress¬ 
ible  or  elastic  in  dilatation.  This  has  been  influenced 
in  part  by  the  fact  that  accurate  measurement  of  visco¬ 
elastic  characteristics  in  dilatation  are  difficult  to 
make,  and  only  now  are  preliminary  results  becoming  avail¬ 
able.  Since  more  precise  measurements  of  dilatational 
response  can  be  expected,  certain  solutions  of  viscoelastic 
stress  analysis  problems  are  presented  to  illustrate  the 
convenience  of  utilizing  measured  relaxation  moduli  in 
dilatation  as  well  as  shear  in  evaluating  stress  distri¬ 
butions.  Examples  for  cylindrical  boundaries  are  given. 

These  solutions  are  also  useful  in  suggesting  experimental 
configurations  from  which  the  relaxation  modulus  in  dilata¬ 
tion  can  beNdeduced  mathematically  while  avoiding  the 
difficulties  Nof  direct  measurement  by  surrounding  a  specimen 
with  fluid  under  pressure.  Metal  casings  can  be  incorporated 
_ 

The  results  communicated  in  this  report  were  obtained 
in  the  course  of  an  investigation  conducted  under  Contract 
Nonr-225(69)  by  Stanford  University  with  the  Office  of 
Naval  Research,  Washington,  D.  C.  Reproduction  in  whole  or 
in  part  is  permitted  for  any  purpose  of  the  United  States 
Government.  This  work  is  to  be  presented  at  the  Fourth 
International  Congress  on  Rheology,  at  Brown  University, 
Providence,  R.  I.,  August  1963- 


into  the  system  to  eliminate  penetration  of  pressurizing 
fluid  into  the  specimen  and  to  permit  more  accurate  strain 
measurement,  and  at  the  same  time  permit  evaluation  of 
the  relaxailon  moduli  of  the  polymer  by  means  of  the 
stress  analysis  theory. 


Introduction 


In  recent  years  rapid  progress  has  been  achieved  in 
the  development  of  methods  of  stress  analysis  for  linear 
isotropic  viscoelastic  bodies  within  the  framework  of 
infinitesimal  strain  and  displacement  theory.  For  such 
bodies,  invariance  of  the  constitutive  equations  under 
rotation  of  axes  implies  that  two  independent  viscoelastic 
operator  relations  are  needed  to  prescribe  the  material 
behavior,  analogous  to  the  two  Independent  elastic  constants 
in  classical  elasticity  theory.  Since,  for  many  visco¬ 
elastic  bodies,  shear  strain  magnitudes  often  dominate 
dilatational  strains,  it  is  appropriate  to  take  as  the  two 
independent  stress-strain  operator  relations  those  govern¬ 
ing  shear  and  dilational  deformation;  for  otherwise  the 
shear  behavior  is  likely  to  dominate  combined  operators 
as  it  does  for  Young's  modulus  for  extension.  The  much 
larger  deformability  in  shear  is  reflected  in  the  literature 
by  the  much  broader  coverage  of  viscoelastic  characteristics 
in  shear  and  the  comparative  paucity  of  dilatational 
response  measurements.  Although  measurements  in  shear 
and  simple  tension  can  formally  be  combined  to  give  the 
separate  shear  and  dilatational  responses,  practically 
the  dilatational  strain  appears  as  the  difference  between 
two  larger  quantities  which  is  therefore  subject  to 
appreciable  error,  and  this  carries  over  to  the  deduced 
bulk  modulus  characteristic  functions.  Direct  bulk  modulus 
measurements  are  difficult  to  carry  out  since  the  specimen 
must  be  subjected  to  hydrostatic  pressure  in  a  cavity,  and 
hence  the  expansion  of  this  cavity,  and  the  compression  of 
the  pressurizing  fluid  must  be  taken  into  account.  Moreover, 
the  fluid  may  penetrate  into  the  specimen,  and  small  trapped 
air  bubbles  can  cause  major  discrepancies  in  measurement. 
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Reliable  response  characteristics  are  needed  for 

satisfactory  stress  analysis,  and  the  lack  of  information 

alluded  to  above  has  hampered  application  of  the  theory, 

12  R 

as  discussed  by  Kolsky  and  Shi  ,  Lee  ,  and  Hunter  .  In 
this  paper  we  re-assess  the  situation,  present  a  method 
of  solution  for  the  stress  distribution  in  a  spinning 
pressurized  circular  cylinder  for  a  material  viscoelastic 
in  dilatation, in  order  to  indicate  that  analytical  methods 
are  available  to  treat  such  problems,  and  show  how  such 
analyses  can  suggest  methods  of  measuring  viscoelastic 
bulk  characteristics. 

Because  of  the  lack  of  information  on  viscoelastic 
characteristics  in  dilatation,  most  treatments  of  visco¬ 
elastic  stress  analysis  problems  have  included  an  assumption 

of  uncertain  validity  to  prescribe  this  property.  Incompressi- 
4  B  6 

bility  is  likely  to  be  satisfactory  if  the  body  is 

of  open  section,  relatively  unconfined  so  that  free  flow 

in  shear  is  permitted.  However,  this  is  likely  to  be  too 

B  7 

drastic  an  assumption  in  many  cases.  Elastic  dilatation-^-’' 
provides  a  more  accurate  representation.  The  latter  intro¬ 
duces  different  operators  in  shear  and  dilatation  and  this 
can  increase  the  complexity  of  certain  types  of  analysis 
as  mentioned  below.  To  avoid  this,  the  assumption  has  been 
made°-’^  that  the  operators  governing  shear  and  dilatational 
response  differ  only  by  a  constant  factor.  This  assumption 
is  equivalent  to  the  condition  that  the  viscoelastic 
operator  corresponding  to  Poisson's  ratio  in  elasticity 
reduces  to  a  constant  as  in  elasticity,  so  that  effectively 
only  one  time  dependent  operator  occurs.  However,  as 
pointed  out  by  Hunter  ,  there  is  little  evidence  to  support 
this  assumption,  although  it  does  provide  a  simple  means 
of  introducing  the  influence  of  viscoelastic  response  in 
dilatation,  albeit  in  an  artificial  manner. 
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The  analytical  convenience  of  the  assumptions  mentioned 
above  is  particularly  significant  for  analyses  using  differ¬ 
ential  operator  stress-strain  relations,  since  the  introduction 
of  viscoelastic  compressibility  would  lead  to  higher  order 
combined  operators  and  hence  call  for  more  Involved  methods 
of  solution.  But  as  .remarked  previously^  ,  the  use  of  the 
Integral  operator  representation  of  viscoelastic  character¬ 
istics  may  not  lead  to  appreciably  more  complicated  analysis 
with  viscoelastic  compressibility,  and  further,  the  technique 
described  in^*^  allows  the  direct  use  of  measured  viscoelastic 
properties,  should  these  be  available  in  the  appropriate 
form.  Hence,  in  the  following  we  shall  incorporate,  after 
suitable  transformation,  the  viscoelastic  dllatational 
characteristics  given  by  Marvin,  Addrich  and  Sack^^  into 
the  stress  analysis  of  a  spinning,  pressurized  cylinder, 

essentially  following  the  same  procedure  as  that  described 
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in  a  previous  analysis  for  an  elastically  compressible 
cylinder. 


Development  of  the  Theory 

We  shall  first  consider  the  theory  for  plane  strain, 

axially  symmetric  viscoelastic  problems.  Apart  from  the 

effect  of  the  different  constitutive  relation  for  dilatation, 

12 

the  analysis  is  similar  to  that  considered  in  ,  In  which, 
the  material  is  considered  to  behave  elastically  under 
hydrostatic  pressure.  Since  plane  strain  conditions  and 
cylindrical  symmetry  are  assumed,  then,  in  conventional 
notation,  the  strains  e^,  e^,  and  the  displacement 
u  are  related  by: 


_  Su  c  -  — 
r  “  3?  ’  ^0  ~  r  ’ 


(1) 


The  stresses  o^,  Oq  and  satisfy  the  equilibrium  equation: 

-  3- 


I 


(2) 


where  p  is  the  mass  density  of  the  viscoelastic  material, 

and  0)  =  cj(t)  is  the  angular  velocity  of  rotation  of  the 
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cylinder.  As  discussed  in  reference  ,  shear  stress  due  to 
angular  acceleration  is  not  considered.  The  stresses  and 
strains  also  satisfy  the  constitutive  equations  of  linear 
viscoelasticity: 


Or  -  Cg  =  2  [  G(t-T)  (e^  -  e^)  dr 


=  2 

r  z 


t 

j  G(t-T)  (£g  -  e^)  dT 


(3) 


(4) 


+  ‘'e  +  °z  =  3 


K(t-T)  ^  (e„  +  Cg  +  €^)  dT  (5) 


This  last  equation  (5)  expresses  the  viscoelastic  behavior 
in  dilatation;  a  material  behaving  elastically  in  dilatation 
is  a  special  case  included  in  (5)  by  making  the  function 
K(t)  a  constant. 

The  material  is  considered  undisturbed  prior  to  zero 

time  so  that  the  lower  limit  -«>  in  the  integrals  in 

(3)-(5)  niay  be  replaced  by  0  .  Then,  just  as  in 
'  1  o 

reference  ,  Eqs.  (l)-(3)  give: 

r  +  pr^‘*)^=  -  2  I  G(t-T)  -  7)  dx  (6) 


I 
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and  hence 
aj.(r,t)  =  f^(t) 

where  Is  an  arbitrary  function  Introduced  by  the 

spatial  Integration.  To  obtain  another  Independent  equation 
relating  and  u,  we  eliminate  Oq  and  from 

Eqs.  (3)-(5)  arid  then  use  (l)  to  give: 

t  t 

3a^(r,t)  =  3f  K(t-T)  +  ^)dT  +  2^  G(t-T)|;^(2|^  -  ^)dT 

'o  ■'o 


1  2  2/,  ^ 
gpr  10  (t) 


G(t-T) 


{r, T)dT 


(7) 


t  2  ^ 

=  [3K{t-T)4G(t-T)]|:^  (ru)dT+3^  G(t-T)|^(^- 

(8) 

In  order  to  leave  the  equation  with  the  u  terms  contained 
exclusively  in  a  complete  derivation,  we  use  Eq.  (6)  to 
eliminate  the  second  integral  in  (8),  then  integrate  with 
respect  to  r  to  give: 


-  ^pr^co^(t)  +  if  [3K(t-T)+G(t-T)  ]|:^eg(r,T)dT 

(9) 


where  fgCt)  is  another  arbitrary  function  Introduced  by 
integration . 

We  therefore  have  two  integral  equations  for  the  four 
unknown  functions  f^(t),  f2{t),  a^(r,t)  and  eg(r,t),  and 
hence  require  two  Independent  boundary  conditions  to 
properly  define  the  problem. 

In  order  to  obtain  the  stress  distributions  It  Is 


•^ic 


convenient  to  obtain  an  equation  not  containing  displace¬ 
ment.  This  is  done  by  operating  on  both  sides  of  Eqs.  (j) 
and  (9)  with  the  linear  integral  operators 

t  t 

[3K(t-T)  +  G(t-T)]|^  dT  and  2f  dx 

r\ 


respectively and  adding  to  give: 


ha  (r,T) 

[3K(t-T)  +  4G(t-T)]-^^^:j^ -  dx 


Sf, (x)  o  f  Sf  (x) 

[3K(t-x)+G(t-x)]— dx+^(  G(t-x)  '  dx 


2  t 

pvj  (6K(t-T)  +  5G(t-T)) 


This  can  be  reduced  to  a  much  more  convenient  form  by 
Introducing  the  auxiliary  function  R(t)  defined  by  the 
integral  equation: 

t 

f  [3K(t-x)  +  4G(t-x)]  dx  =  2G(t)  (11) 

•^6" 

Then  rearranging  Eq.  (lO)  so  as  to  include  all  K  terms 

only  in  the  form  3K  +  4G  ,  and  using  (ll),  linearity, 

13 

and  the  associative  property  of  convolutions  , 
cfjo(r,t)  is  given  by  the  simpler  equation: 


-6- 


t)  =  f^(t)  -  I  pr'^u^(t) 


+  |(  [-^ 


fi(T)  *£^pA]  ax 


Which  may  be  rewritten  as: 


a  (r,t)  =  P  (t)  +  ^^-P  (t)r^  • 
j-  r 

where  P^(t)  =  R(t-T)  ^ 


o/  ,  Sf2(T) 
Ro(t)=-?l  R(t-T)  — dT 


dr 


(13) 

(I4a) 

(14b) 


P3(t)  [4w2(t)  -  sf  R(t-T)  dt]  (I4c) 


using  (13)  and  (2),  Oq  is  given  by: 

Pp(t)  P  P 

ag(r,t)  =  P,(t) - 2 —  +  [3P3(t)  +  pu'^(t)]r  (15) 


Por  the  non-rotating  cylinder  (0=0)  Eqs .  (l3)and  (15) 
are  immediately  recognized  as  the  Lame  solution,  which  is 
independent  of  material  properties  if  the  two  independent 
boundary  conditions  specify  the  surface  tractions.  Por 
the  more  general  rotating  case  (l4c)  gives  P2(t)  and 
knowledge  of  the  internal  and  external  pressures  then 
immediately  determines  Pj^(t)  and  P2(t)  and  hence  the 


(12) 
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stress  distributions.  If,  however,  one  or  both  surfaces 

has  a  bonded  elastic  casing,  the  stress  solution  -is  not 

so  simple.  At  such  a  surface,  the  boundary  condition 

14 

relates  the  normal  stress  to  the  displacement  according 
to; 


(b,t)  =  -  B  eg(b,,t) 


(16) 


where  b  is  the  boundary  radius,  and 


1 

B 


1  + 


yhil  -  vl) 


%{l  +  v^) 


(17) 


where  Ej^,  and  h  are  the  Young's  modulus,  Poisson's 
ratio  and  thickness  respectively  of  the  reinforcing  elastic 
casing. 

If  only  one  surface  has  a  bonded  elastic  casing,  and 

the  other  has  prescribed  normal  pressure,  the  method  of 

12 

solution  is  precisely  that  described  in  reference  with 
the  R-function  now  defined  by  (ll).  The  method  Involves 
using  the  pressure  boundary  condition 

^^.(a^t)  =  - 

and  the  casing  boundary  condition  (l6)  together  with 
(7)  for  r  =  b,  and  (12)  for  r  =  a,b  to  eliminate  f2(t) 
and  obtain  two  interrelated  convolution  integral  equations 
for  a^(b,t)  and  f^(t),  which  can  be  simultaneously 
solved  numerically  by  finite  difference  techniques. 
cf^(r,t)  and  aQ(r,t)  are  then  immediately  given  by 
Eqs.  (13)  and  (I5). 

If  both  boundary  conditions  are  of  the  casing  type  (16), 
then  both  can  be  applied  to  each  of  Eqs.  (7)  and  (9) 
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permitting  elimination  of  fj^(t)  and  f2(t)  and  hence 
again  obtain  two  interrelated  convolution  integral  equations, 
but  this  time  for  o^{a,t)  and  a^(b,t)  (or  the  equivalent 
boundary  strains).  Either  f^(t)  or  f2(t)  are  then  given 
by  simple  quadratures,  and  substitution  in  Eqs.  (13)  and 
(15)  determines  a^(r,t)  and  aQ(r,t). 

So  far  we  have  only  considered  the  determination  of  the 
stress  distributions.  Obtaining  the  displacement  field 
would  seem  to  involve  only  a  simple  choice  between,  solving 
the  Integral  equations  (7)  or  (9)  for  eg(r,t)  and  hence 
u(r,t).  Of  these  two  equations,  (7)  seems  to  be  the 
preferable  since  it  may  be  rewritten  in  the  quadrature  form; 


eg(r,t)  =  if  J{t-T)  [f^(T)  -  -  CT^(r,T)]  dx 


where  J(t)  is  the  creep  compliance,  either  a  measured 
material  characteristic  or  determined  from  G(t)  as  by 
Hopkins  and  Hamming  However,  in  order  to  achieve 
satisfactory  numerical' accuracy,  Eq.  (9)  is  in  fact  to 
be  preferred  and  when  rewritten  in  the  form: 


[3K(t-T)  +  G(t-T)]  [b^eg(b,T)  -  r^eg(r,T)]  dr 


=  0  (b,t)  -  ra{r,t)]  +  4(b^-r^)pw^(t )  (20) 


it  may  be  solved  as  a  Volterra  Integral  equation  for  e^fr,!). 
Another  equation  for  eg(r,t)  ,  not  containing  G(t) 
directly,  may  be  obtained  by  eliminating  thfe  G  term  from 
Eq.  (7)  and  (9),  and  f2(t)  through  the  substitution 
r  =  b,  to  give: 
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t 


6K(t-T)  |:p  [b^eg(b,T)  -  r^e0(r,T)]  dt  =  b^[4a^(b,t)-f^(t)  ] 

-  r^[4a^(r,t)-fj_(t)] 

(21) 


This  equation  is  more  suitable  for  the  determination  of 
displacements  for  elastically  compressible  materials  since 
it  then  reduces  to  a  simple  algebraic  expression  for  eg(r,t). 

Particular  Example 

The  particular  example  to  be  considered  in  this  section 
is  that  of  a  hollow  cylinder  of  polyisobutylene  encased  in 
an  elastic  shell  and  subjected  to  arbitrary  internal  pressure. 
Although  in  the  numerical  example  presented  later  the 
cylinder  is  non-rotating  and  non-ablating,  the  theory 
developed  in  the  previous  section,  and  the  method  of 
numerical  evaluation  used,  apply  for  the  general  case. 

The  choice  of  material  and  loading  conditions  was 
governed  by  the  availability  in  the  literature  of  measured 
physical  characteristics  in  shear  and  dilatation.  Marvin, 
Aldrich  and  Sack^^  give  complex  shear  and  bulk  moduli  for 
polyisobutylene  over  an  equivalent  frequency  range  of 
10^  to  10^  cycles/sec.  when  the  measurements  are  reduced 
to  25°C.  Although  these  cannot  be  substituted  directly 
into  the  theory  presented  here,  which  is  based  on 
relaxation  moduli,  in  the  absence  of  measurements  of  the 
latter,  we  decided  to  base  the  solution  on  the  relaxation 
moduli  deduced  approximately  from  the  complex  modulus 
measurements,  since  our  objective  was  to  present  an 
example  illustrative  of  the  method  which  might  serve  to 
encourage  more  precise  measurements  of  viscoelastic 
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dllatational  characteristics  to  permit  application  of  this 
approach.  We  therefore  used  the  empirical  formulae 
proposed  by  Ninomiya  and  Perry^^  for  evaluating  approximate 
relaxation  moduli: 


G(t)  =  G'(u)  -  0.40  G"  (0.40io)  +  0.014  G"  (lOo) 
K(t)  =  K'(u))  -  0.40  K"  (0.40a))  +  0.0l4  K"  (lOw) 


li)  = 


1/t 


=  1/t 


(22a) 


(22b) 


Pig.  1  shows  the  resulting  moduli  plotted  logarithmically, 

the  experimental  results  having  been  smoothed  by  differencing. 

Because  of  the  equivalent  high  frequency  at  the  reduced 

temperature  of  25°C,  direct  evaluation  of  the  relaxation 

moduli  by  (22)  would  cover  a  time  range  in  micro-seconds, 

which  would  be  unsuitable  for  the  quasi-static  stress 

analysis  problem  under  consideration.  We  therefore  consider 

pressurizing  to  occur  at  a  lower  temperature  which  intro- 

4 

duces  a  factor  of  10  in  the  time  scale  on  the  basis  of 

17 

the  thermo-rheologically  simple  hypothesis  .  This 
corresponds  to  a  temperature  of  about  ^4  which 

quasi-static  analysis  would  be  satisfactory  for  loading 
times  in  the  range  of  the  relaxation  data.  The  time  scale 
in  Fig.  1,  and  the  example  to  be  presented,  correspond 
to  this  temperature,  although  the  solution  can  be  inter¬ 
preted  to  correspond  to  different  loading  times  at  different 
temperatures,  as  long  as  the  temperature  is  not  considered 
raised  to  the  extent  that  inertia  forces  due  to  the 
deformation  become  significant.  It  is  seen  from  Pig.  1 
that  much  more  marked  relaxation  occurs  in  shear  than  in 
dilatation,  but  that  the  latter  occurs  sooner, 
r  R(t)  was  obtained  from  the  integral  equation  (11 ) 

by  finite  sum  numerical  evaluation  as  described  in 
reference^^.  Fig.  2  shows  the  variation  of  R(t) 
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computed  using  several  step  lengths.  It  Is  evident  that 
practical  convergence  is  achieved  for  A(log  t)  =  O.O625. 
For  comparison  the  corresponding  curves  for  elastically 
compressible  bodies  are  shown,  assuming  the  constant 
bulk  modulus  to  coincide  with  the  initial  modulus,  K(0), 
and  long-time  modulus,  K(<»),  respectively.  Comparison  of 
the  influence  of  the  time  step  A(  log  t)  with  that  for  the 
elastically  compressible  case^*^  reveals  much  slower  con¬ 
vergence  for  viscoelastic  dilatation.  Whereas  in  the 
elastically  compressible  case  R(t)  is  monofronically 
decreasing.  Pig.  2  shows  this  to  be  no  longer  so.  That 
this  behavior  is  permissible  is  clear  from  a  physical 
interpretation  of  the  meaning  of  the  function. 

R(t)  defined  by  (11 )  is  a  generalization  of  the 

7 

R  function  Introduced  by  Muki  and  Sternberg  .  It  has 
a  physical  interpretation  in  terms  of  one-dimensional 
strain  with  displacement  parallel  to  the  axis  and 

Independent  of  the  Xg,  x^  coordinates.  The  principal 
strains  are  then  (e^,  0,  0)  and  stresses  (a^,  c^), 

and  the  linear  viscoelastic  laws  in  shear  and  dilatation 
give : 


c^l(t)  -  ^^(t)  =(  2  G(t-T)  37^  dT 


(23a) 


(J^(t)  +  2c72(t)=( 


3K(t-T) 


(23b) 


Elimination  of  <3^  yields: 


/  oe, 

3a^(t)  =  [3K(t-T)  +  4G(t-T)]  ^  dt 
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Comparison  with  (ll)  shows  R(t)  to  be  the  strain  response 

e^(t)  with  this  one-dimensional  constraint  for  a  stress 

variation  2G(t)/3,  and  hence  €2(t)  can  increase  with 

time  as  does  creep  response. 

In  the  previous  section  we  showed  that  for  a  cavity 

of  radius  a  subjected  to  pressure  p(t),  with  an  elastic 

casing  at  r  =  b^  a  pair  of  Integral  equations  for 

(j^(b,t)  and  f^(t)  can  be  deduced  for  (7)  and  (12). 

Transformed  into  a  convenient  form  for  numerical  analysis 

12 

they  reduce  to  the  equations  given  in  reference  for 
the  elastically  compressible  case,  modified  to  accommodate 
the  minor  change  in  the  definition  of  R: 

t 

[1-|  G(0)]  a^(b,t)  +  |f  G(t-T)dT  =  f^(t)-ipr^u)2(t) 

4  (25a) 

t 

[l-^R(0)Jf^  (t)  +  -^  (  R(t-T)dT 

■'o 

=  [b^a^(b,t)+a4(t)+(b^-a4Fo(t)]  (25b) 

b"^-a 

12 

As  described  in  reference  these  can  be  solved  conveniently 
by  finite  sum  approximation,  in  a  manner  similar  to  that 
for  single  integral  equations^^.  Programming  was  carried 
out  in  the  computer  language  BALGOL,  for  processing  on  the 
IBM  7090  machine  at  the  Stanford  University  Computing 
Center.  These  facilities  were  in  part  made  available 
under  the  National  Science  Foundation  Grant  NSP-GP  9^8. 

In  order  to  have  an  appreciable  influence  of  visco¬ 
elasticity  in  dilatation,  Pig.  1  shows  that  a  loading 
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time  of  the  order  10  ^  secs,  is  needed,  and  so  the  pressure 
variation  in  the  cavity  was  taken  to  be: 

p(t)  -  p_{l  -  (26) 

Since  the  stiffness  of  the  external  elastic  casing 

is  likely  to  have  a  significant  influence  on  the  stress 
12 

distributions  ,  the  solution  was  carried  out  for  a 
"soft"  casing,  B  =  1.434  GCo),  which  corresponds  to 
an  aluminum  casing  with  i  ~150^  a  "medium"  casing, 

B  =  3'441  G(0),  which  corresponds  to  a  steel  casing  with 
~8o,  and  for  a  rigid  casing. 

Table  I  shows  the  dimensionless  compressive  stress 
values  (a/p^)  for  various  times  and  boundary  positions 
for  the  medium  casing,  compared  with  the  corresponding 
problems  with  elastic  compressibility  and  K  =  K(0)  and 
K(“)  respectively.  Because  of  the  rapid  initial  variation, 
equal  steps  of  log  t  were  taken  for  the  finite  sum 
approximation  for  t  <  0,01,  and  equal  steps  of  t  there¬ 
after.  The  steps  were  halved  until  a  difference  of  less 
than  in  resulted,  this  occurring  for  the  visco¬ 

elastic  case  with  A(log  t)  =  0.125  and  At  =  O.OO5.  The 
results  in  Table  I  show  that  for  the  longer  times  the 
elastically  compressible  solution  with  K  =  K(«>)  provides 
a  good  approximation  for  both  stress  components,  as  would 
be  expected  from  the  form  of  the  relaxation  modulus  in 
dilatation  shown  in  Pig.  1.  Viscoelasticity  in  dilatation 
influences  the  component  appreciably  at  the  shortest 

time,  and  particularly  the  value  of  0^(0,!).  At  the 
later  times,  the  dominating  relaxation  in  shear  leads  to 
almost  hydrostatic  stress  conditions  with  the  two  stress 
components  almost  equal. 

For  the  softer  casing  larger  displacements  occur,  with 


consequent  larger  stress  differences,  and  for  the  same 
time  steps  the  effect  of  halving  the  step  length  produced 
a  change  of  about  2%  in  showing  somewhat  larger, 

though  acceptable,  truncation  errors. 

It  was  pointed  out  at  the  end  of  the  previous  section 
that  displacements  could  be  evaluated  from  the  Inversion  of 
(7)  in  the  form  (19)  which  depend  on  the  viscoelastic  shear 
operator,  or  by  (20)  which  contains  an  operator  dominated 
by  the  dilatational  viscoelastic  kernel  K(t)  (since  as 
shown  in  Pig.  1,  K(t)»G(t)).  For  an  encased  cylinder, 

after  the  initial  rapid  growth  of  displacement,  continued 
expansion  is  inhibited  by  the  casing  and  the  dilatational 
material  characteristics,  and  shear  restraint  virtually 
disappears  as  the  stress-differences  relax.  Thus  when 
using  (19)  the  constraining  effect  must  arise  since 

1  p  p 

fl(T)  -  -^pr  u)  -  o^(r,T)  ~  0.  In  the  example  considered 

here,  u  =  0,  and  the  computation  of  the  displacement 
then  depends  on  the  difference  of  two  nearly  equal  com¬ 
puted  quantities.  This  leads  to  considerable  numerical 
error,  and  as  shown  by  the  points  plotted  in  Pig.  3  for 
0^(3, t)  =  for  the  medium  casing,  B  =  Cr(0), 

satisfactory  convergence  was  not  achieved  even  for  time 
increments  one  quarter  those  used  for  the  stress  calculations. 
Use  of  (20),  in  which  the  constraint  on  u  is  introduced 
in  a  natural  way  through  the  dilatational  operator,  no 
such  difficulty  arose  and  the  change  in  u  for  halving 
the  time  step  for  the  courser  increments  used  for  the 
stress  calculations  was  less  than  1^. 

Pig.  3  shows  the  circumferential  strain  at  the  cavity 
surface,  eg(a,t)  =  u(a,t)/a,  thus  being  a  measure  of 
displacement,  for  the  two  elastic  casings  mentioned  above, 
and  for  a  rigid  casing.  Values  calculated  for  elastic 
compressibility,  for  K(0)  and  K('»),  are  also  shown. 
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In  Pig.  3  Eq  Is  presented  in  dimensionless  form 

£q  =  G(0)  —  which  accounts  for  the  magnitude  of  order 

unity.  It  is  evident  that  the  stiffness  of  the  casing 

has  a  major  effect  on  the  displacements^  and  that  with 

increasing  time  the  displacement  for  the  body  exhibiting 

viscoelasticity  in  dilatation  approaches  the  elastically 

compressible  case  with  K  =  K(<»).  This  would  be  expected 

in  view  of  the  K(t)  variation  shown  in  Pig.  1,  and 

the  influence  of  dilatational  constraint  in  this  problem. 

The  approach  however  is  slow,  and  in  fact  the  limit  would 

only  be  achieved  when  a  hydrostatic  stress  condition  has 

developed.  Table  I  shows  that  this  has  not  been  achieved 

at  t  =  0.1  since  (a^  -  Og)  is  still  appreciable,  and 

a  can  be  expected  to  differ  by  an  amount  of  the  same 

z 

order.  Only  when  the  maximum  shear  stresses,  l.e.  the 
stress  differences,  have  relaxed  to  zero,  will  elastic 
dilatational  constraint  based  on  K(oo)  determine  the 
displacement  correctly. 

Por  the  rigid  casing  example,  incompressible  material 
behavior  would  give  u  h  0,  and  we  see  that  the  elastically 
compressible  case  based  on  K(0)  is  closer  to  this  result 
than  to  the  dilatational  viscoelastic  and  K('»)  cases. 

Por  K(o),  the  Poisson's  Ratio  of  instantaneous  response 
is  0.43,  compared  with  0.27  for  K(<»).  It  is  thus  clear 
why  the  former,  with  growth  towards  0.5  as  the  shear 
modulus  relaxes,  is  close  to  the  incompressible  solution. 

The  results  discussed  above  indicate  that  in  this 
problem  dilatational  viscoelasticity  has  overall  only  a 
moderate  Influence  on  the  stress  and  displacement  distri¬ 
butions.  If  such  a  configuration  were  envisaged  for 
measuring  dilatational  properties,  therefore ,  facilities 
for  accurate  measurements  of  surface  displacement  would  be 
needed.  This  consideration  is  assessed  in  the  following  section. 
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Application  for  Measuring  Dllatatlonal  Characteristics 


The  analyses  given  above  to  determine  stresses  and 
displacements  for  materials  exhibiting  viscoelasticity 
In  dilatation  as  well  as  In  shear,  can  be  used  In  reverse 
to  determine  the  relaxation  moduli  In  shear  and  dilatation 
from  displacement  measurements.  However,  the  comparative 
stiffness  In  dilatation  calls  for  extremely  accurate 
displacement  measurements  to  achieve  satisfactory 
dllatatlonal  modulus  values.  Measurement  of  surface 
displacements  or  strains  for  viscoelastic  bodies  is 
hampered  by  the  texture  of  many  polymeric  specimens 
which  have  a  geometrically  less  smoothly  defined  surface 
than  for  structural  metals,  and  which  are  often  too  soft 
compared  with  resistance  strain  guage  elements,  or  modified 
by  glue  used  for  attaching  them,  to  permit  sufficient 
accuracy  of  measurement.  If  a  cylindrical  specimen  can 
be  reinforced,  both  Internally  and  externally,  by  elastic 
metal  casings,  accurate  displacement  measurements  could 
be  made  on  the  metal  surfaces,  following  Internal  pressure 
or  spin,  and  the  analysis  given  above  can  be  transformed 
for  the  determination  of  G(t)  and  K(t). 

For  example.  If  we  consider  the  Internal  and 
external  shell  constants  (17)  to  be  and  B,  respectively, 
and  consider  expansion  due  to  spin  only,  equation  (20)  can 
be  expressed  In  the  form 


L(t-T)  dT 


Y(t) 


where  L(t)  =  3K(t)  +  G(t) 

X(t)  =  b  u(b,t)  -  a  u(a,t) 


(27) 


Y(t)  =  -  3[b  u(b,t)-aBg  u  (a  ,  t )  ]/2+3  (b^-a^  )p  a)^(t  )/8 
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and  equation  (7)  in  the  form: 


t 

2I  G(t-T)  =  V(t)  (28) 

'0 

where  U(t)  =  u(b,t)/b  -  u(a,t)/a 

V(t)  =  -  B^u(b,t)/b  -  BgU(a,t)/a  -  pa)2(t )  (b2-a^)/2 

For  measured  u(b,t)  and  u(a,t)  ,  (27)  and  (28)  can  be 

solved  numerically  for  L(t)  and  G(t)  just  as  (ll)  was 
solved  for  R(t)^^.  K(t)  then  follows  from: 

K(t)  =  [L(t)  -  G(t)]/3  (29) 

Measurements  of  the  displacement  on  steel  or  aluminum 
rings,  for  example,  may  prove  sufficiently  accurate  to 
permit  satisfactory  evaluation  of  K(t). 

Similar  approaches  can  be  adopted  for  other  loading 
situations  which  may  be  more  convenient  experimentally 
than  the  above  suggestion  based  solely  on  centrifugal 
Inertial  loading.  This  was  chosen  for  illustrative 
purposes  to  permit  direct  application  of  equations  already 
developed  in— this  paper,  Eqs.  (16)  and  (17)  can  be  readily 
modified  to  include  internal  pressure  on  the  inner  ring, 
as  well  as  support  pressure  at  the  metal-plastic  interface, 
so  that  a  similar  integral  equation  for  the  evaluation  of  K(t) 
follows  from  displacement  measurements  following  appli¬ 
cation  of  a  known  Internal  pressure  variation.  Such  a 
configuration  would  have  the  advantage  of  maintaining 
pressure  at  the  metal-plastic  Interfaces  and  hence  avoid 
adhesion  problems.  Alternatively,  for  known  internal 
pressure  without  an  inner  ring,  measurements  of  the 


displacement  of  the  outer  casing  could  be  used  to  determine 
K(t‘J,  if  G(t)  were  measured  independently ^  for  example 
by  measurements  of  torsion.  Since  the  method  of  evaluation 
is  numerical  solution  of  integral  equations,  direct  com¬ 
putation  from  experimental  readings  of  displacement  or 
pressure  is  available. 

A  similar  but  simpxer  analysis  can  be  carried  out  for 
the  spherically  symmetrical  problem,  simpler  since  only 
two  independent  principal  stress  components  need  be 
considered  in  the  theory  ( o^,  Oq)  compared  with  three 

in  the  cylindrical  case.  However,  the  construction  of 
spherical  specimens  with  well  fitting  elastic  casings  is 
likely  to  present  extremely  difficult  problems  of  fabri¬ 
cation,  and  hence  it  seems  unlikely  that  this  approach 
will  lead  to  applications.  The  analysis  is  therefore  not 
presented . 
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Table  I  Compressive  Stress  Values 


,  ^-1000 
1-e  J 


t  -> 

0.001 

0.01 

0.1 

K(0) 

.489 

.949 

a^(b,t)  K(t) 

.460 

.920 

K(co) 

.452 

.918 

K(0) 

.260 

.867 

.978 

ag(b,t)  K(t) 

.207 

.792 

.962 

K(oo) 

.163 

.788 

.964 

a^(a,t)  all  K 

.632 

1.000 

l.OOC 

K(0) 

.118 

.817 

.96s 

CTg(a,t)  K(t) 

.043 

.712 

.948 

K(oo) 

.018 

.673 

.95c 
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Log  (G  or  K  dynes /cm^) 


in  Dilatation  (K)  and  Shear  (G) 


Log,p  (t  secs) 


Fig.  2  Characteristic  Functions  R(t)  for  Polyisobutylene 


Fig.  3  Variations  of  the  Internal  Boundary  Displacements  for 
different  Compressibilities  and  elastic  casings 


